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SUMMARY 


A  unidirectional  fiberglass  epoxy  composite  specimen  is  modelled 
as  a  homogeneous  transversely  isotropic  continuum  plate-like  medium. 
Acousto-ultrasonic  non-contact  input-output  characterization  by  trac¬ 
ing  the  P  waves  in  the  continuiom  is  studied  theoretically  with  a 
transmitting  and  a  receiving  transducer  located  on  the  same  face  of 
the  plate.  The  isotropic  plane  of  the  equivalent  continutam  plate 
model  lies  in  the  midplane  of  the  plate  and  is  parallel  to  the  top  and 
the  bottom  faces  of  the  plate. 

The  single  reflection  problem  for  an  incident  P  wave  at  a  stress - 
free  plane  boundary  in  a  semi- inf inite  transversely  isotropic  medium 
whose  isotropic  plane  is  parallel  to  the  plane  boundary  is  analyzed 
first.  It  is  found  that  an  obliquely  incident  P  wave  results  in  a 
reflected  P  wave  and  a  reflected  SV  wave.  One  quadrant  of  each  of  two 
sheets  of  the  slowness  surfaces  of  the  reflected  P  and  SV  waves  is 
plotted.  The  angle  of  reflection  of  the  reflected  P  wave  is  equal  to 
the  angle  of  incidence  of  the  incident  P  wave.  However,  the  angle  of 
reflection  of  the  reflected  SV  wave  is  smaller  than  the  angle  of  inci¬ 
dent  of  the  incident  P  wave.  Consequently,  no  critical  angle  phenome¬ 
non  occurs.  The  amplitude  ratios  of  the  reflected  P  and  SV  waves  to 
the  incident  P  wave  as  a  function  of  the  angle  of  incidence  are 
plotted.  The  balance  in  energy  flux  normal  to  the  plane  boundary  is 
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checked. 


The  delay  time  for  propagation  between  the  transmitting  and  the 
receiving  transducers  is  computed  as  if  the  P  waves  were  propagating 
in  an  infinite  half  space.  It  is  found  that  the  directional  depen¬ 
dence  of  the  phase  velocity  of  the  P  waves  propagating  in  the  trans¬ 
versely  isotropic  medium  has  a  significant  effect  on  the  delay  time, 
as  opposed  to  the  directional  independence  of  the  phase  velocity  of 
the  P  waves  propagating  in  an  isotropic  medium. 


The  displacements  associated  with  the  P  waves  in  the  plate  and 
which  may  be  detected  by  the  non-contact  receiving  transducer  are 
approximated  by  an  asymptotic  solution  for  an  Infinite  transversely 
isotropic  medium  subjected  to  a  harmonic  point  load.  The  polar  dia¬ 
grams  for  the  directivity  functions  are  plotted  at  frequencies  of 
0.75,  1.50  and  2.25  MHz. 

The  study  enhances  the  quantitative  understanding  of  acousto- 
ultrasonic  nondestructive  evaluation  (NDE)  parameters  such  as  the 
stress  wave  factor  (SWF)  and  wave  propagation  in  fiber  reinforced  com¬ 
posites  or  any  other  materials  which  can  be  modelled  as  transversely 
isotropic  media. 
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INTRODUCTION 


Fiber  reinforced  composite  materials  are  attractive  materials  for 
aerospace  applications  because  of  their  high  specific  mechanical  prop¬ 
erties.  It  has  been  shown  that  many  composites,  such  as  fiberglass 
epoxy  composites  or  fiber  reinforced  ceramics,  as  shown  in  Fig.  1,  may 
be  modelled  as  homogeneous  transversely  isotropic  continue  [1].  In 
this  work,  acousto -ultrasonic  (AU)  non-contact  input-output  character¬ 
ization  of  a  homogeneous  transversely  isotropic  elastic  plate  is 
investigated  by  tracing  P  waves. 

First,  the  single  reflection  problem  of  an  incident  P  wave  at  a 
stress-free  plane  boundary  in  a  semi- inf inite  transversely  isotropic 
medium,  whose  isotropic  plane  is  parallel  to  the  plane  boundary  is 
considered.  At  such  boundaries,  the  conditions  for  the  existence  of 
wave  mode  conversion,  the  angle  of  reflection  of  the  reflected  wave, 
and  the  amplitude  ratios  of  the  reflected  wave  to  the  incident  wave 
are  derived. 

Second,  the  P  wave  input- output  relations  are  derived  when  mul¬ 
tiple  reflections  occur  at  the  top  and  bottom  faces  of  the  plate.  The 
delay  time  between  input  and  output  versus  the  distance  separating  the 
transmitting  and  receiving  transducers  is  analyzed.  The  directivity 
functions  of  the  stresses  associated  with  the  P  waves  are  computed. 
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And,  the  output  displacement  at  the  non- contact  receiving  transducer 
is  approximated  by  an  asymptotic  Solution. 

This  investigation  should  enhance  the  quantitative  understanding 
of  AU  NDE  parameter  such  as  the  stress  wave  factor.  It  also  provides 
the  potential  for  assisting  in  the  development  of  better  NDE  schemes 
utilizing  the  SWF  as  well  as  other  AU  parameters . 
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STNGLE  REFLECTION  PROBLEM  AT  STRESS-FREE  PLANE  BOUNDARY  IN  SEMI- 
INFINITE  TRANSVERSELY  ISOTROPIC  MEDIUM  WHOSE  ISOTROPIC  PLANE  IS  PARAL¬ 
LEL  TO  PLANE  BOUNDARY  FOR  INCIDENT  P  WAVE 

1 .  Reflected  P  and  SV  Waves 

For  a  homogeneous  linearly  elastic  transversely  isotropic  contin¬ 
uum,  the  niomber  of  independent  elastic  constants  is  five  [1].  Define  a 
coordinate  system  (x,  y,  z)  for  a  semi-infinite  transversely  isotropic 
medium  whose  isotropic  plane  is  parallel  to  the  plane  boundary  where 
the  reflection  occurs  as  follows:  the  plane  boundary  contains  the  x 
and  y  axes,  and  the  z  axis  is  the  zonal  axis  of  the  medium,  which  is 
in  the  direction  parallel  to  the  fiber  direction  shown  in  Fig.  1.  See 
Fig.  2.  The  generalized  Hooke's  law  is  written,  relative  to  the  (x, 
y,  z)  coordinate  system,  as  [1] 

’•xx  -  Ciiu,x  +  Ci2V,y  +  Ci3W,z 
’■yy  “  +  ^ll'''^>y  C]^3W,z 

Tzz  ”  Ci3U,x  -I-  Ci3V,y  +  C33W,z  (1) 

Txz  ”  C44(u,z  +  w,x) 

ryz  =  C44(v,z  +  w,y) 

^xy  “  C56(u,y  -i-  v,x) 

where  Tj-g  (r,  s  -  x,  y  and  z)  are  the  normal  (r  -  s)  and  shear  (r  s) 
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stresses  with  respect  to  the  coordinate  system  (x,  y,  z) ;  u,  v  and  w 
are  the  displacement  components  of  a  point  in  the  medium  along  the  x, 
y  and  z  axes,  respectively;  denotes  partial  differentiation  with 

respect  to  the  variable  which  follows;  and  C]L2’  ^13 »  ^33  ^44 

are  the  five  independent  elastic  constants  where  ==  1/2  (Cn  -  ^12^ 


Let  a  plane  progressive  wave  be  represented  as  [2] 

(u,  V,  w)  =  A  (P^,  Py,  P^;)  expliwCS^x  +  Syy  +  S^z  -  t)  }  (2) 

where  S^,  Sy  and  are  the  components  of  the  slowness  vector,  which 
is  in  the  same  direction  as  the  normal  to  the  wavefront  and  whose  mag 
nitude  is  equal  to  the  reciprocal  of  the  magnitude  of  the  phase  velo¬ 
city  [1],  along  the  x,y  and  z  axes,  respectively;  P^,  Py  and  P^  are 
the  components  of  a  unit  vector  of  particle  displacement  along  the  x, 
y  and  z  axes,  respectively;  A  is  the  amplitude  of  particle  displace¬ 
ment;  t  denotes  time  and  u)  denotes  radian  frequency. 


It  follows  from  Eqs .  (1)  and  (2)  that  the  stresses 
sented  as 

=  iwA[C;LlSxPx  ^I2^y^y  Ci2^S^?^]exip[ic^(S^yL  +  Syy  + 

Tyy  =  iwA[Cl2SxPx  +  Cl^SyPy  +  C  ^  3  S  ^  P  Z  ]  6  X?  {  i  CO  (  S  +  Syy  + 

=  iwA[C]_3SxPx  +  ClsSyPy  +  C33S2PZ ]  exp { icoCS^x  +  Syy  + 

Txz  =  i«A[C44S2Px  +  C44SxPz]exp{ico(SxX  +  Syy  +  S^z  -  t)  } 

ry2  =  icoA[C44S2Py  +  C44SyP2]exp{iw(SxX  +  Syy  +  S^z  -  t) } 

r^y  “  icoA[C56SyPx  +  CggS^PylexpliwCS^x  +  Syy  +  S^z  -  t) ) 


can  be  repre 

S^z  -  t)) 

S^z  -  t) ) 

S^z  -  t)} 

(3) 
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The  stress  boundary  conditions  on  the  stress -free  plane  boundary 


require  that  [2] 


^xz 

^yz 

’’zz 


(I) 

(I) 

(I) 


+  r 


xz 


+  T 


yz 


+  r 


zz 


(R)  -  0 
(R)  -  0 
(R)  =  0 


(4) 


where  (r=x,  y  and  z)  represents  stresses  on  the  plane  boundary 

associated  with  the  incident  P  wave,  and  (r  =  x,  y  and  z) 

r0pres6nts  stresses  on  the  plane  boundary  associated  with  the 
reflected  waves. 

In  order  to  satisfy  Eq.  (4),  it  is  required  that  the  frequency, 
u)y  of  the  reflected  waves  be  equal  to  that  of  the  incident  wave  and 
that  [ 2 ] 


(R) 


Sy(I)  - 


(5) 


As  a  result  of  Eqs .  (5),  the  slowness  vectors  of  the  incident 

and  reflected  waves  lie  in  a  plane  called  the  plane  of  incidence. 

This  analysis  can  be  simplified  by  assuming  that  the  plane  of  inci¬ 
dence  is  the  X  =  0  plane;  that  is,  the  slowness  vectors  of  the  inci¬ 
dent  and  reflected  waves  are  in  the  x  =  0  plane,  as  shown  in  Fig.  2. 
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Then  it  follows  from  Eqs .  (5)  that 


8^(1)  =  =  0. 


(6) 


It  has  been  shown  [3]  that  P  waves  or  SV  waves  travelling  in  a 
plane  containing  the  zonal  axis,  z  axis,  of  a  transversely  isotropic 
medium  are  quasi -longitudinal  and  quasi -  transverse ,  respectively;  that 
is,  the  components  of  the  unit  vector  of  particle  displacement  for 

either  the  P  or  SV  waves  along  both  the  y  and  z  axes,  P^  and  P  do 

y  z 

not  vanish;  whereas  the  components  along  the  x  axis,  P^,  do  vanish. 
Therefore,  it  follows  from  Eqs.  (3)  and  (6)  that  stresses  associated 
with  P  and  SV  waves  are 


^xy  ’’xz 

=  0 

ryz  ^  ^  \ 

;  rxx  ^  0  ; 

’■yy  ^  0  ; 

^ZZ  ^  ^ 

(7) 


It  has  also  been  shown  that  an  SH  wave  travelling  in  a  plane 
containing  the  zonal  axis,  z  axis,  of  a  transversely  isotropic  medium 
possesses  a  transverse  displacement  only;  that  is,  for  the  coordinates 
in  Fig.  2  (P^^,  P^,  P^)  =  (1,0,0)  [3].  Therefore,  it  follows  from 
Eqs,  (3)  and  (6)  that  the  stresses  associated  with  SH  waves  are 


T 

T 


XX 

xz 


=  Tyy 
^  0  ; 


’’zz  ^yz 


Txy  ^  0. 


0 


(8) 


8 


Assume  that  a  P  wave  is  incident  on  the  plane  boundary.  It  fol¬ 


lows  from  Eqs.  (4)  and  (7)  that 


r  (I) 
^yz 


^  0; 


zz 


(I)  ^  0 


xz 


(I)  =  0. 


(9) 


As  a  result  of  Eq.(9),  it  is  immediately  known  from  Eq.  (4) 
that  is  equal  to  zero.  This  means  that  no  SH  wave  will  be 

reflected  back  into  the  medium  because  a  reflected  wave  of  the  SH 
type  would  result  in  nonzero  values  of  the  stress  So,  Eqs.  (4) 

reduce  to: 


^  (I)  +  r  -  0 

^yz  ^  ^yz  ^ 

^  (I)  +  r  (R)  -  0. 

’^ZZ  ’zz 


(10) 


Since  either  a  reflected  P  wave  or  a  reflected  SV  wave  results 
in  nonzero  values  of  the  stresses,  and  it  is  therefore 

concluded  from  Eqs.  (7)  and  (10)  that  both  a  P  wave  and  an  SV  wave  may 
be  reflected  back  into  the  medium. 


2 .  Slowness  Surface  for  P  and  SV  Waves 

The  equations  of  motion  relative  to  the  coordinate  system 
(x,y,z) 


9 


are  [1] 


XX  ,  X 

’’xy.y 

^XZ  ,  z 

“  pu, tt 

xy,x 

+ 

’’yy.y 

+ 

^yz,z 

-  pv,tt 

(11) 

xz  ,x 

’■yz.y 

+ 

^zz ,  z 

=  pw,tt 

where  the  body  forces  are  identically  zero  for  the  homogeneous  solu 
tion. 


It  follows  from  Eqs .  (1),  (2)  and  (11)  that  the  following  equa- 
tions  of  motion  are  obtained: 

[CiiS^^  +  C66Sy2  +  C44S22  -  p]F^  +  (C12  +  C66)SxSyPy 

+  (Ci3  +  C44)SxSyP2  =  0 

<^12  +  CgglSxSyPx  +  +  CiiSy2  +  C44S22  -  p]Py  (12) 

+  (Ci3  +  C44)SyS2P2;  =  0 

(Ci3  +  C44)SxSyPx  +  (Ci3  +  C44)SyS^Py  +  [044(5^2  +  s^2) 

+  -P]Pz  =  0 

The  condition  for  the  existence  of  the  plane  wave  solution  is 
expressed  by  setting  the  determinant  of  the  matrix  of  the  coefficients 
of  Pj^,  Py  and  P^  in  Eq.  (12)  equal  to  zero  [1]: 
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[CllSx^  *  C66Sy2  +  .^3 

(C12  *  C64)S^Sy 

(Ci3  +  C^4)SxSz 

<^^12  *  ''66^®x®y 

CC66Sx2  *  C33S/  ♦  -/,] 

(C33  +  C^4)SySj, 

(Ci3  +  C44)SySj 

[C44(S^2  ^  5^2)  ^  c33S^2  .^3 

(13) 

By  expanding  Eq.  (13),  three  sheets  of  slovmess  surface  are 
obtained.  The  slowness  surface  for  an  P  wave  is  given  in  [3]  as 


Cl  1 +C 


ll'^^44  2  2  ‘^44''’‘^33  2 

- )  (Sx+S„)  +  ( - ) 

2  2 


2  2 

+  1/2  {[(Cii-C44)(Sx+Sy  ) 


+  (C33  -C44)S^]2 

2  2  2 

+  4(Sjj+Sy)S2[  (C2^2^-C44)  (C33-C44) 


-(013+044)^]}^/^  “  P’  (1^) 

Similarly,  the  slowness  surface  for  an  SV  wave  is  given  in  [3] 


C11+C44 


)(S"+Sv)  +  (- 


C44+C33  2 


-  1/2  {[(Cii-C44)(S^+S^)  +  (C33-C44)s2]2 

+  Hsl+sl)sl  [(Cii-C44)(C33-C44) 


-((Ci3+C44)2])V2  =  P 


n 


In  the  present  study,  the  slowness  vectors  of  the  incident  P 
wave  and  the  reflected  P  and  SV  waves  are  confined  in  the  y-z  plane, 
as  shown  in  Fig.  2.  Take  numerical  values  of  elastic  constants  and 
density  given  in  [1]  for  the  unidirectional  fiberglass  epoxy  composite 
as  follows:  C]_]_  =  10.581  x  10^  N/m^ ,  C|^3  =  4.679  x  10^  N/m^,  C33  - 
40.741  X  10^  N/m^,  C44  =  4.422  x  10^  N/m^ ,  and  p  =  1850  kg/m^ .  One 
quadrant  of  the  intersection  of  the  slowness  surface  of  a  P  wave 
travelling  in  the  unidirectional  fiberglass  epoxy  composite  with  the 
plane  x  =  0  and  one  quadrant  of  the  intersection  of  the  slowness  sur¬ 
face  of  an  SV  wave  travelling  in  the  unidirectional  fiberglass  epoxy 
composite  with  the  plane  x  -  0  are  obtained  by  substituting  the  numer¬ 
ical  values  of  the  elastic  constants  and  the  density  given  above  into 
Eqs.  (14)  and  (15),  respectively,  and  are  shown  in  Fig.  3. 

3 ,  Angle  of  Reflection 

It  follows  from  Eqs.  (5)  and  (6)  that  the  y-component  of  the 
slowness  vector  of  an  incident  P  wave  is  equal  to  the  y-component  of 
the  reflected  P  wave  as  well  as  the  y-component  of  the  reflected  SV 
wave.  Accordingly,  the  relation  between  the  y-component  of  the  slow¬ 
ness  vector  of  an  incident  P  wave  and  the  y-components  of  the  slowness 
vectors  of  the  reflected  P  and  SV  waves  is  given  as 

Syd)  -  Sy<‘')  -  SyfSV)  .  b  (16) 
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where  represents  the  y-component  of  the  slowness  vector  of  an 

incident  P  wave;  Sy^^^  represents  the  y-component  of  the  slowness 
vectors  of  the  reflected  P  wave;  Sy^^^)  represents  the  y-component  of 
the  slowness  vector  of  the  reflected  SV  waves;  and  b  is  a  common  con¬ 
stant  . 


It  follows  from  Eqs .  (6),  (14)  and  (16)  that  the  relation 
between  the  z- component  of  the  slowness  vector  of  an  incident  P  wave 
and  that  of  the  reflected  P  wave  is 

(17) 

The  minus  sign  is  due  to  the  fact  that  the  slowness  vector  of  an 
incident  P  wave  points  out  of  the  medium,  whereas  the  slowness  vector 
of  the  reflected  P  wave  points  into  the  medium,  as  shown  in  Fig.  2. 
Consequently,  the  value  of  the  z-component  of  the  slowness  vector  of 
an  incident  P  wave,  S2(^\  is  negative,  whereas  that  of  the  reflected 
P  wave,  S2(^\  is  positive. 

The  angle  of  reflection  is  defined  as  the  angle  between  the 
slowness  vector  of  a  reflected  wave,  either  type  P  or  SV ,  and  the 
normal  to  the  plane  boundary  where  the  reflection  occurs.  Similarly, 
the  angle  of  incidence  is  defined  as  the  angle  between  the  slowness 
vector  of  an  incident  P  wave  and  the  normal  to  the  plane  boundary,  as 
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shown  in  Fig.  2.  Therefore,  the  angle  of  reflection  of  a  reflected  P 
wave  is 

=  tan-l(Sy(P)/S2(^))  (18) 

and  the  angle  of  incidence  of  an  incident  P  wave  is 

$1  =  tan-l(Sy(I)/-S2^^)).  (19) 

It  follows  from  Eqs .  (16),  (17),  (18)  and  (19)  that  the  angle 

of  incidence  of  an  incident  P  wave  is  equal  to  the  angle  of  reflection 
of  the  reflected  P  wave,  as  shown  in  Fig.  2. 

However,  the  angle  of  reflection  of  the  reflected  SV  wave  is  not 
equal  to  the  angle  of  incidence  of  the  incident  P  wave.  For  a  given 
value  of  b  in  Eq.  (16),  two  values  of  the  z-component  of  the  slowness 
vector,  S^,  of  equal  magnitude  but  opposite  sign  of  an  SV  wave  travel¬ 
ling  in  the  plane  x  =  0  in  the  transversely  isotropic  medium  can  be 
obtained  from  the  slowness  surface  for  the  SV  wave  from  Eqs.  (6)  and 
(15) .  The  positive  z-component  of  the  slowness  vector  corresponds  to 
the  reflected  SV  wave  and  is  denoted  as  Similarly,  for  a 

given  value  of  b,  there  exists  a  positive  z-component  of  the  slowness 
vector,  S^^^^  corresponding  to  the  reflected  P  wave.  Accordingly, 
for  any  given  value  of  b  in  Eq.  (16),  there  exist  a  positive  z-compo- 
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nent  of  the  slowness  vector  for  the  reflected  P  wave,  and  a 

positive  z-component  of  the  slowness  vector  for  the  reflected  SV  wave, 

In  fact,  the  z- components  of  the  slowness  vectors  of  the 
reflected  P  and  SV  waves,  S2  (P)  and  S2(SV)_  for  an  incident  P  wave 
which  determines  the  value  of  b  in  Eq.  (16)  and  which  travels  in  the 
plane  x  =  0  in  the  transversely  isotropic  medium,  can  be  obtained  from 
the  lengths  of  the  perpendicular  lines  between  the  abscissa  represent¬ 
ing  the  value  of  b  and  the  interesections  with  two  sheets  of  slowness 
surface  for  the  reflected  P  wave  and  the  reflected  SV  wave,  respec¬ 
tively,  as  shown  in  Fig.  3.  It  is  apparent  from  Fig.  3  that  the 
z-component  of  the  slowness  vector  of  the  reflected  P  wave,  is 

less  than  that  of  the  reflected  SV  wave,  ,  for  an  incident  P 

wave  travelling  in  the  plane  x  -  0  in  the  unidirectional  fiberglass 
epoxy  composite  under  consideration  [1].  Consequently,  the  angle  of 
reflection  of  the  reflected  SV  wave  defined  similarly  to  Eq.  (18) 

as 

*SV  -  «n-l  (Sy(SV)/s^<SV),,  (20) 

is  smaller  than  the  angle  of  reflection  of  the  reflected  P  wave,  from 
Eqs.  (16),  (18)  and  (20).  The  angle  of  reflection  of  the  reflected  SV 
wave  as  a  function  of  the  angle  of  incidence  of  an  incident  P  wave 
travelling  in  the  unidirectional  fiberglass  epoxy  composite  [1]  is 
shown  in  Fig.  4. 
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The  critical  angle  is  defined  as  the  angle  of  incidence  of  an 
incident  P  wave  at  which  the  slowness  vector  of  the  reflected  SV  wave 
(or  P  wave)  becomes  tangent  to  the  plane  boundary  where  the  reflection 
occurs.  For  the  case  of  an  incident  P  wave  travelling  in  the  unidirec¬ 
tional  fiberglass  epoxy  composite  in  Fig.  1,  it  has  been  shown  that 
the  angle  of  reflection  of  the  reflected  SV  wave  is  smaller  than  that 
of  the  reflected  P  wave.  Thus,  the  critical  angle  phenomenon  does  not 
occur  for  the  unidirectional  fiberglass  epoxy  composite  under  consid¬ 
eration.  This  is  due  to  the  fact  that  when  the  angle  of  incidence  of 
an  incident  P  wave  reaches  90°,  the  angle  of  reflection  of  the 
reflected  SV  wave  is  still  less  than  90°. 


Amplitude  Ratios  of  Reflected  Waves  to  Incident  Wave 


It  has  been  shown  that  when  a  P  wave  travelling  in  a  semi- 
infinite  transversely  isotropic  medium  is  incident  on  a  plane  bound- 
ary ,  a  P  wave  and  an  SV  wave  will  be  reflected.  The  boundary  condi¬ 
tions  on  the  stresses  of  a  P  wave  travelling  in  the  plane  x  =  0  and 
incident  on  the  plane  boundary  can  be  obtained  from  Eqs ,  (4)  and  (7) 
as 


+  ,„(SV)  +  ,„(P)  .  0 


zz 


zz 


(21) 
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where  and  represent  the  shear  and  normal  stresses 

associated  with  the  incident  P  wave;  represent 

the  shear  and  normal  stresses  associated  with  the  reflected  SV  wave; 

and  r  and  represents  the  shear  and  normal  stresses 

yz  z  z  * 

associated  with  the  reflected  P  wave. 

The  shear  stress,  associated  with  a  P  wave  of  unit 

amplitude  travelling  in  the  plane  x  =  0  and  incident  on  the  plane 
boundary  at  the  origin,  as  shown  in  Fig.  2,  can  be  obtained  from 
Eqs .  (3)  and  (6)  as 

-  i»(C55S^(OPy<I)  +  C55Sy<I)p^(I))  (22) 

where  S  and  are  the  components  of  the  slowness  vector  of  the 

y  ^ 

incident  P  wave  along  the  y  and  z  axes,  respectively;  Py^^^  and  ^ 

are  the  components  of  the  unit  vector  of  particle  displacement  of  the 
incident  P  wave  along  the  y  and  z  axes,  respectively;  C55  is  an  elas¬ 
tic  constant;  and  cj  is  the  frequency.  Similarly,  the  normal  stress, 
associated  with  the  incident  P  wave  of  unit  amplitude  can  be 
expressed,  from  Eqs,  (3)  and  (6),  as 

.  MCi3Sy<I)Py<l)  +  C33S^<»P^<I))  (23) 

where  C^3  and  C33  are  elastic  constants . 
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The  shear  stress,  and  the  normal  stress, 

associated  with  the  reflected  P  wave  on  the  plane  boundary  at  the 
gin  are,  from  Eqs .  (3)  and  (6), 


r  (P)  -  WP)(C55S,(P)P  (P)  +  C55S  (P)P^(P)) 


la,A<P)(Ci3Sy(P)Py<P)  +  C33Sj(P)p,(P)). 


where  and  are  the  components  of  the  slowness  vector  of  the 

reflected  P  wave  along  the  y  and  z  axes,  respectively;  Py^^^  and 
P^^^^  are  the  components  of  the  unit  vector  of  particle  displacement 
of  the  reflected  P  wave  along  the  y  and  z  axes,  respectively;  and 
a(P)  is  the  amplitude  of  the  reflected  P  wave. 


The  shear  stress,  ^nd  the  normal  stress, 

associated  with  the  reflected  SV  wave  on  the  plane  boundary  at  the 
origin,  from  Eqs,  (3)  and  (6), 


+  C55S  (26) 


r^ySV)  ,  i^A(SV)(Ci3Sy(SV)py(SV)  ^  C33Sy  SV)  (SV)  )  (27) 

where  Sy^SV)  are  the  components  of  the  slowness  vector  of 
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the  reflected  SV  wave  along  the  y  and  z  axes,  respectively  ;  Py(SV)  and 
p^(SV)  are  the  components  of  the  unit  vector  of  particle  displacement 
of  the  reflected  SV  wave  along  the  y  and  z  axes,  respectively;  and 
^(SV)  amplitude  of  the  reflected  SV  wave. 

Upon  substitution  of  Eqs.  (22)  through  (27)  into  Eq.  (21),  the 
boundary  conditions  on  the  stresses  for  an  incident  P  wave  of  unit 
amplitude  travelling  in  the  plane  x  =  0  can  be  rewritten  as 


(Ci3Sy(nPy<n  t  C33S,(I)p^(I)) 

+  A(f)(Ci3Sy<f)Py<'')  +  C33Sj(f)Pj('’>)  (28) 

+  A<SV)(Ci3Sy(SV)Py(SV) 

+  C33Sj(SV)p^(SV))  .  0 

The  components  of  a  unit  vector  of  particle  displacement  along 
the  y  and  z  axes  of  a  P  wave  travelling  in  any  plane  containing  the 
zonal  axis  of  a  transversely  isotropic  medium  are  given  in  [3],  and 
can  be  expressed,  when  applied  to  the  present  case  of  the  reflected  P 
wave  travelling  in  the  plane  x  -  0,  as  [3] 
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Py(P)  =  [Hp  -  (C33-C44)S2^P)^/ 

{[Hp-(C33-C44)S2(P)Y 

+  [(Ci3+C44)S 

(29) 

P^^P)  =  (Ci3+C44)Sy(P)s2(PV 

([Hp-(C33-C44)S^(P)Y 

+  [(Cl3+C44)Sy(P)s2^P)]^^^^ 

where  C]_3,  C33  and  are  elastic  constants;  and  Hp  is  defined  as 

Hp  =  {(Cii-C44)Sy(P)^  +  (C33-C44)S2^P)^ 

+  {[(Cii-C44)Sy(P)^ 

+  (C33-C44)S2(P)^^ 

2  2 

-  4Sy<P)  [(Cii-C44)(C33-C44) 

-  (Ci3+C44)^]^)^/^)/2- 


Similarly,  the  components  of  a  unit  vector  of  particle  displacement 
along  the  y  and  z  axes  for  the  reflected  SV  waves  travelling  in  the 
plane  x  ==  0  can  be  expressed  as 
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I 


Py<^''>  -  [Hsv  -  (C33-C44)Sj(SV)2j/ 
l[Hsv  ■  «:33-C44)Sz<®''>  1 

+  ((Ci3+C44)Sy(S'')s,(SV)]^^/^  (30) 

P,(SV)  .  (Ci3+C44)Sy<SV)s^(SV)/ 

IlHsv  -  «:33-':44)Sz<®''>  1^ 

+  [<Cl3+C44)Sy<=''>S,<SV)jS^/^ 

where  is  defined  as 

Hsv  =  {(Cii-C44)Sy<SV)2  ^ 

-  ([(Cu-C44)Sy(SV)2  +  (C33-C44)S,(SV)Y 

-  4Sy(SV)2  g^(sv)2  [(Cii.C44)(C33-C44) 

2  2  1/2 

-  (C13+C44)  ]  )  )/2. 

The  amplitude  ratios  of  the  reflected  P  and  SV  waves,  and 

a(SV),  for  a  P  wave  of  unit  amplitude  obliquely  incident  on  the  plane 
boundary  can  be  determined  from  Eq.  (28).  For  a  given  incident  P 
wave,  the  values  of  the  components  of  the  slowness  vector  along  the  y 
and  z  axes,  ^y^^^  values  of  the  components  of  the 

unit  vector  of  particle  displacement  along  the  y  and  z  axes,  a^d 
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are  defined  as  part  of  the  the  specification  of  the  incident  P 
wave.  The  values  of  the  components  of  the  slowness  vector  of  the 
reflected  P  wave  along  the  y  and  z  axes,  and  are 

determined  from  Eqs .  (16)  and  (17),  respectively.  The  values  of  the 
components  of  the  unit  vector  of  particle  displacement  along  the  y  and 
z  axes  for  the  reflected  P  wave,  Py^^^  and  are  obtained  by 

substituting  Sy^^)  and  into  Eq .  (29).  The  values  of  the 

coefficients  for  in  Eq.  (28)  are  thus  obtained.  Since  the  y-com- 

ponent  of  the  slowness  vector  of  the  reflected  SV  wave,  is 

equal  to  that  of  the  incident  P  wave,  Eq .  (16),  the  value  of  the 
z-component  of  the  reflected  SV  wave,  is  obtained  by  substi¬ 

tuting  Sy^^)  into  the  slowness  surface  for  the  SV  wave,  Eq.  (15).  On 
substitution  of  the  values  of  the  and  Sy^^^)  into  Eq .  (30),  the 

values  of  the  unit  vector  of  particle  displacement  along  the  y  and  z 
axes  for  the  reflected  SV  wave,  Py^^^^  and  are  determined. 

The  values  of  the  coefficients  for  in  Eq .  (28)  are  thus 

obtained.  The  values  of  the  amplitudes  of  the  reflected  P  and  SV 
waves,  A^^)  and  A^^^\  are  then  obtained  by  solving  Eq.  (28)  with  the 
thus  determined  values  of  the  coefficients  for  A^^^  and  A^^^)  in  Eq. 
(28)  for  a  given  incident  P  wave.  By  varying  the  angle  of  incidence 
of  the  incident  P  wave  of  unit  amplitude  and  by  repeating  the  proce¬ 
dures  described  above,  the  amplitude  ratios  of  the  reflected  P  and  SV 
waves  to  the  incident  P  wave  are  obtained  as  functions  of  the  angle  of 
incidence.  For  an  incident  P  wave  in  the  unidirectional  fiberglass 
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epoxy  composite  shown  in  Fig.  1,  the  amplitude  ratios  of  the  reflected 
P  and  SV  waves  to  the  incident  P  wave  versus  the  angle  of  incidence 
are  shown  in  Fig.  5. 

5 .  Balance  in  Energy  Flux  Normal  to  Plane  Boundary 

The  balance  in  energy  flux  normal  to  the  plane  boundary  z  =  0, 
as  shown  in  Fig.  2,  must  be  satisfied  and  is  expressed  as  [2] 

F^(I)  +  F^(P)  +  F^CSV)  _  0  (31) 

where  and  are  the  z-components  of  the  energy 

fluxes  of  the  incident  P  wave,  the  reflected  P  wave  and  the  reflected 
SV  wave,  respectively. 

The  z- component  of  the  energy  flux  of  an  incident  P  wave  of  unit 
amplitude  travelling  in  the  plane  x  —  0  is  [4] 

F2(I)  =  «2(C44Py(I)^  +  Ci3Py(I)P2(^)Sy(I) 

+  C44Py(I)P2<^)Sy(I)  +  C33P2^^>^S2^^^)  (32) 

where  Py^^^  and  P2^^^  the  components  of  the  unit  vector  of 

particle  displacement  of  the  incident  P  wave  along  the  y  and  z  axes, 
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respectively;  and  are  the  components  of  the  slowness 

vector  of  the  incident  P  wave  along, the  y  and  z  axes,  respectively;  w 
is  the  radian  frequency;  and  €2^3,  C33  and  C44  are  elastic  constants. 

Similarly,  the  z- components  of  the  energy  fluxes  of  the 
reflected  P  wave  and  the  reflected  SV  wave  are  [4] 

F^(P)  .  a(P)2  ■o2(C44Pj,<P>^  S,(P)  +  Ci3Py(P)p^(P)Sy<P) 

+  C44Py<P)p2<P)Sy<P)  +  C33P2<P)^Sj(P))  (33) 

and 


F^(SV)  .  a(SV)2^2(c^^P  (SV)2s^(SV)  +  Ci3P„(SV)p^(SV)s  (SV) 


+  C44P  (SV)p^(SV)s  (SV)  ^  C33P,(SV)2s^(SV)) 


(34) 


Since  the  frequency  term  is  common  to  Eqs .  (32),  (33)  and 
(34),  the  balance  in  energy  flux  normal  to  the  plane  boundary, 

Eq.  (31),  is  not  affected  by  assuming  the  value  of  the  frequency  w  to 
be  equal  to  unity.  Accordingly,  subsequent  calculations  of  the  values 
of  the  z- components  of  the  energy  fluxes  of  the  incident  P  wave  and 
the  reflected  P  and  SV  waves,  and  are  done  by 

assuming  the  radian  frequency  w  to  be  equal  to  one. 
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The  values  of  the  z- components  of  the  energy  fluxes  in 
Eqs.  (32),  (33)  and  (34)  are  obtained  similarly  to  the  calculations 
of  the  amplitude  ratios  of  the  reflected  P  and  SV  waves  to  the  inci¬ 
dent  P  wave.  For  a  given  incident  P  wave  of  unit  amplitude,  the 
y-components  of  the  slowness  vectors  and  of  the  unit  vectors  of  par¬ 
ticle  displacement,  ?y^^\  ^y^^^ 

the  z- components  of  the  slowness  vectors  and  of  the  unit  vectors  of 
particle  displacement, 

of  the  incident  P  wave  and  of  the  reflected  P  and  SV  waves  are  deter¬ 
mined  first.  Then,  combining  the  amplitudes  of  the  reflected  P  and  SV 
waves,  and  with  the  values  of  ,  Sy^^^ ,  Sy^^^^ , 


values  of  the  z- components  of  the  energy  fluxes  of  the  incident  P  wave 
and  the  reflected  P  and  SV  waves,  and  are  thus 

obtained  from  Eqs.  (32),  (33)  and  (34).  For  a  P  wave  of  unit  ampli¬ 
tude  travelling  in  the  plane  x  -  0  in  the  unidirectional  fiberglass 
epoxy  composite  shown  in  Fig.  1  incident  on  a  plane  boundary,  the 
z- components  of  the  energy  fluxes  of  the  reflected  P  and  SV  waves  and 
the  energy  flux  of  the  incident  P  wave  are  shown  in  Fig.  6  with  the 
value  of  the  frequency  w  in  Eqs.  (32),  (33)  and  (34)  equal  to  one.  The 
balance  in  energy  flux  normal  to  the  plane  boundary,  Eq.  (31),  is  also 
checked  and  is  shown  in  Fig,  6. 
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ACOUSTO-ULTRASONIC  NON-CONTACT  INPUT-OUTPUT  CHARACTERIZATION  OF  UNIDI¬ 


RECTIONAL  FIBERGLASS  EPOXY  COMPOSITE  PLATE 

It  has  been  shown  [1]  that  the  unidirectional  fiber  composite 
shown  in  Fig.  1  may  be  modelled  as  a  homogeneous  transversely  iso¬ 
tropic  continuum.  For  the  axes  shown  in  Fig.  1,  the  isotropic  plane 
of  its  equivalent  continuum  lie  in  the  midplane  of  the  plate  [1].  A 
cartesian  coordinate  system  (x,  y,  z)  is  chosen  so  that  the  x-y  plane 
is  the  isotropic  plane;  as  a  result,  the  upper  and  the  lower  surface 
are  at  z  -  h/2  and  z  -  “h/2,  respectively,  where  h  is  the  plate  thick¬ 
ness.  The  properties  of  the  equivalent  continuum  model  of  the  unidi¬ 
rectional  fiberglass  epoxy  composite  plate  to  be  considered  are  [1] 

h  -  0.1  ra 

Cii  -  10.581  X  10^  N/m2 

C]^3  -  4.679  X  10^  N/m2 

C33  -  40.741  X  10^  N/m2  (35) 

C44  -  4.422  X  10^  N/m2 

Cgg  -  3.243  X  10^  N/m^ 

p  -  1850  kg/m^ 

A  transmitting  and  a  receiving  transducer  are  located  on  the 
same  face  of  the  unidirectional  fiberglass  epoxy  composite  plate  spe¬ 
cimen  without  direct  contact,  as  shown  in  Fig.  7.  The  unidirectional 
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fiberglass  epoxy  composite  plate  specimen  shown  in  Fig.  7  is  consid 
ered  as  a  plate  of  thickness  h  and  of  infinite  planar  (x-y)  extent. 

The  input  electrical  voltage  to  the  transmitting  transducer  is  V^Ct) 
and  the  output  electrical  voltage  from  the  receiving  transducer  is 
Vo(t)  where  t  represents  time.  The  transmitting  transducer  converts 
an  input  electrical  voltage  into  a  stress,  whereas  the  receiving 
transducer  converts  a  displacement  associated  with  stress  waves  trav¬ 
elling  in  the  plate  into  an  output  voltage.  In  the  following  analy¬ 
sis,  only  the  P  waves  are  traced.  The  P  waves  which  are  generated  by 
the  transmitting  transducer  located  above  point  0  experience  multiple 
reflections  at  each  face  of  the  plate,  and  then  reach  the  receiving 
transducer  located  above  point  M,  as  shown  in  Fig.  8.  Since  the  iso¬ 
tropic  plane  lies  in  the  midplane  and  is  parallel  to  both  the  top  and 
the  bottom  faces  where  the  multiple  reflections  occur,  the  angle  of 
reflection  of  the  reflected  P  wave  is  equal  to  the  angle  of  incidence 
of  an  incident  P  wave  for  each  reflection  at  each  face  of  the  plate . 
Accordingly,  the  P  wave  travelling  from  the  input  0  to  the  output  M 
may  be  considered  as  a  wave  propagating  in  a  semi- infinite  trans¬ 
versely  isotropic  medium,  and  travelling  to  point  M'  as  if  there  were 
no  bottom  face,  as  shown  in  Fig.  8. 

1.  Delay  Time  and  Phase  Velocity 

Let  the  input  0  and  the  output  M  lie  in  the  y-z  plane.  Assume 
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the  number  of  reflections  at  the  bottom  face  experienced  by  the  P 
wave  in  travelling  from  the  input  0  to  the  output  M  is  n,  as  shown  in 
Fig,  8.  With  respect  to  the  z  axis,  the  angle  of  incidence  of  the  P 
wave  at  each  face  of  the  plate  is  0,  and  the  total  distance  travelled 
by  the  wave  is  R^.  From  the  geometry  in  Fig.  8, 

e  =  tan‘^(i/2nh)  (36) 

where  i  is  the  separation  distance  between  the  input  0  and  the  output 
M,  and 

-  Jl/sin$  (37) 

The  time  delay  t^  for  the  wave  to  reach  the  receiving  transducer  is 

=  Rn/Cl(0  (38) 

where  C]_(^)  is  the  directionally  dependent  phase  velocity  of  the  P 
wave.  The  phase  velocity  Cj_  of  a  P  wave  in  the  unidirectional  fiber¬ 
glass  epoxy  composite  shown  in  Fig.  1  is  given  as  [1] 

Ci(^)  =  [(C44+C2^1  sin^^  +  C33  cos'^B  +  jE)/2p]'^/^  (39) 

where 

E  ”  nCii-C44)  sin^^  +  (044-033)003^6 ] ^ 

+  '^(C]^3+044)^  sin^6  cos^d  ; 
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Cii,  C]^3,  C33,  C44  and  p  are  given  by  Eq.  (35). 


The  delay  time  is  then  computed  when  the  number  of  reflections  n 
at  the  bottom  face  of  the  plate  is  equal  to  10,  100,  300  or  500.  The 
numerical  results  are  shown  in  Fig.  9  where  the  delay  time  t^  is 
plotted  as  the  ordinate,  and  the  dimensionless  separation  i/h  is 
plotted  as  the  abscissa,  for  values  of  zero  to  300.  The  phase  velo¬ 
city  as  a  function  of  the  angle  of  incidence  9  is  also  shown  in 
Fig.  10. 

2 .  Displacements  Detected  by  Receiving  Transducer 

The  displacements  detected  by  the  non-contact  receiving  trans¬ 
ducer  above  point  M,  radiated  by  the  non- contact  transmitting  trans¬ 
ducer,  are  assumed  to  be  equivalent  to  the  displacements  at  point  M' 
associated  with  the  P  wave  propagating  in  a  semi- inf inite  transversely 
isotropic  medium  as  if  there  were  no  bottom  boundary  (except  for  the 
cumulative  effect  of  the  reflection  coefficient),  as  shown  in  Fig.  8. 
The  displacements  at  point  M'  are  approximated  by  the  far- field  asymp¬ 
totic  solution  for  large  of  an  infinite  transversely  isotropic 
medium  subjected  to  a  harmonic  point  load. 

Consider  an  infinite  transversely  isotropic  medium  in  which  the 
z  axis  of  a  rectangular  cartesian  system  0(x,  y,  z)  is  the  zonal  axis 
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of  the  medium  and  the  x-y  plane  coincides  with  the  isotropic  plane,  as 
shown  in  Fig.  11.  The  equations  of  motion  including  the  body  forces 
are  [5] 

^xx,x  ^xy,y  "*■  ^xz ,  z  +  =  pu,tt  (40) 

^xy,x  ^yy , y  ^yz,z  ^  ^  pv,tt  (41) 

^xz,x  ^  ^yz,y  ^zz,z  ==  pw,tt  (42) 

where  (r,  s  =  x,  y  and  z)  are  the  normal  (r  =  s)  and  shear  (r  s) 

stresses  with  respect  to  the  chosen  coordinate  system  0(x,  y,  z) ;  u, 

V  and  w  are  displacement  components  of  a  point  in  the  medium  along  the 
X,  y  and  z  axes,  respectively;  X,  Y  and  Z  are  the  components  of  the 
body  force  along  the  x,  y  and  z  axes  respectively;  p  is  the  density;  t 
is  time;  and  "  denotes  partial  differentiation  with  respect  to  the 
variable  which  follows. 


where  F  and  A  are  given  by  [5]. 
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r  -  w,z;  A  -  u,x  +  v.y. 


By  differentiating  Eq.  (40)  with  respect  to  x  and  Eq.  (41)  with 
respect  to  y,  we  find  upon  addition  of  the  resulting  equations  and 
using  the  appropriate  stress-strain  relations,  Eq.  (1),  that 

Cl t+Ca A  ^44 

A,tt  =  ( — - )  (r,xx  +  r.yy)  + -y 

+  — (A, XX  +  A.yy)  +  X,x  +  Y,y  (44) 

P 

For  a  harmonic  point  load  at  the  origin,  the  body  forces  may  be 
taken  of  the  form  [5] 

X  =  Xo5(x)6(y)5(z)e-i‘^t 

Y  -  Yo5(x)5(y)«(z)e-i‘^^  (^5) 

Z  -  Zo5(x)5(y)«(z)e-i"'' 

where  5(r)  (r  -  x,  y  and  z)  is  the  Dirac  delta  function;  and  Xq ,  Yq 
and  Zq  are  the  magnitudes  of  the  respective  point  body  forces. 

Express  F  and  A  as  threefold  Fourier  integrals  [5]  as  follows. 
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00 


r(x,y,z.t)  =  r  (Sx.Sy.S^.t) 


exp{iw(SxX+Syy+S2Z-t) IdS^dSydS^ 


A(x,y,z,t)  =  ACSj^.Sy.S^.t) 


expliwCS^x+Syy+S^z-t) JdS^dSydS^ 


where 


r(Sx,  Sy ,  ,  t)  =  r(x,y,z,t) 


exp{iw(-SxX-Syy-S2Z+  t))dxdydz. 


A(Sx,  Sy ,  Sz  ,  t)  =  l/8?r^  A(x,y,z,t) 


exp{iw(-Sj^x-Syy-S2Z+  t))dxdydz. 


Similarly,  X,x,  Y,y  and  Z,z  can  be  expressed  as  threefold 
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Fourier  integrals  [5] 


X,x  -  iwSx  X  exp{iw(SjjX  +  Syy  +  S^z  -  t)}dSxdSydS2  (^7) 


Y,y  -  iwSy  Y  exp{iw(SxX  +  Syy  +  S^z  -  t))dSxdSydS2  (48) 


Z,z  -  iwS^  Z  exp{iw(SxX  +  Syy  +  S^z  -  t))dSxdSydS2  (49) 


where 


X  -  l/8ir^  Xo5(x)5(y)5(z)e- 


exp{iw(-SjjX  -  Syy  -  S^z  +  tjdS^^dSydS^ 


Xq/Stt- 


and,  similarly, 


_  _  ‘'O 

Y - -  ;  Z - - 


Substitution  of  Eqs.  (45)  through  (49)  into  Eqs .  (43)  and  (44) 
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In  fact,  Sy,  S^)  =  0  represents  two  sheets  of  the  slowness 

surface,  one  for  a  P  wave  and  one  for  an  SV  wave  [5].  As  a  result  of 
Eqs .  (50)  and  (51),  F  and  A  in  Eq.  (46)  can  be  written  as 
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r  (x,y,z,t)  - 


f  Sz  (is  Xo+iS  Yq) 

JJ  [  87r3wH(Sx.S  S^)  ^ 

Gll/P  (Sx^+Sy^)  C33/P  5^2  .  1  ^ 

8»r^wH(Sx,  Sy ,  S2) 

exp  { iwCSj^x+Syy+S^z- 1)  IdS^dSydS^ 


A(x,y,z, t) 


87r^a)H(S^,  Sy ,  S^) 


(iS^X^+iSyY^) 


-  (CI3+C44)//)  (Sjj2+Sy2)  _ 


8jr-^wH(Sjj^,  Sy ,  $2) 


exp  {iwCS^x  +  Syy  +  S2Z  -  t) )dSxdSydS2  (54) 

The  asymptotic  solution  at  a  large  distance  from  the  point  load  is 
obtained  by  applying  the  theory  of  residues,  the  method  of  stationary 
phase,  and  the  radiation  condition  [5]  as 


^n*^  -(C13+C44)  ^*2,,^*  ^  ^  , 

r(X,y,Z,t)  - -  -  S2  (iSx  Xg+lSy  Yq) 

27rR  p 


*  (£y  (S*%  sJ2)  sf  -  1)1S*  z„ 


exp{ia)(S^x  +  Syy  +  S^z  -  t) } 


[  C44 

2^rR  1  ^ 

*2  *2. 

(S/+  Sy^)  +- 

C33  *2 

P  2 

1  -  (Cl3+‘^44) 

P 

-  (S*^+  Sy^) 

*  \ 

*  * 

* 

exp{iw(S^x  +  Syy  + 

S^z  -  t)} 

^x  o  y  O'' 
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where  R  is  the  distance  from  the  origin  0  where  the  point  loads  ar 

applied  to  the  location  of  interest  Q  in  the  medium,  as  shown  in 
•k  -k  -k 

Fig.  11;  (Sjj,  Sy,  $2)  are  points  on  the  slowness  surface  for  a  P  wave 
where  the  normal  is  parallel  to  the  OQ  direction;  is  the  amplitude 
coefficient  and  is  given  by 


where  E  is  the  sum  with  respect  to  cyclic  permutation  of  S^,  Sy  and 

k  k  'k 

S^,  and  is  evaluated  at  points  (Sj^,  Sy,  S^)  on  the  slowness  for  a  P 
wave  where  the  normal  is  parallel  to  the  OQ  direction. 


The  displacement  components  along  the  x,  y  and  z  axes,  u,  v  and 
w,  due  to  a  P  wave  can  be  obtained  by  direct  integration  of  the  defi¬ 
nitions  of  r  and  A  in  Eq.  (43),  and  are  given  as  [6] 


1 

w - r 

iwS^ 
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Substitution  of  Eqs.  (55)  and  (56)  into  Eq.  (58)  gives  the 
asymptotic  solutions  at  a  large  distance  of  the  displacement  compo¬ 
nents  along  the  x,  y  and  z  axes,  u,  v  and  w,  as  follows: 


-iAn 


27rR(S*^+  Sy^) 


.C44  ,  *2^  *2.  *2  Tw.c:*  Y  .  Y  ^ 

(S^  +  Sy  )  +  ^  -l](iS^  Xq  +  iSy  Yq) 


“(Ci3+Caa)  *2  ■*2  * 

+ - £2 — ^  +  Sy^)  iS^  Zq 


ic  "k  if 

exp{iw(S5^x  +  Syy  +  S^z 


■  '  —  -  |[£^(S*2  +  S*^)  S*^-  l](iS*  Xq  +  iSy  Yq) 

27rR(S*^+  Sy^)  ^  ^  ^ 


-(C13+C4A)  *2  *2,.  * 

- ■---—  (S/4-  Sy^)lS2  Z, 


exp{iw(S^x  +  Syy  +  S^z  -  t) } 


-(Ci3  +  Caa)  *  *  *  V 

- — —  $2  (Sx  Xq  +  Sy  Yq) 


.  sf.  1, 


"k  k  k 

exp{ia)(SjjX  +  Syjr  +  S^z  -  t) 
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The  slowness  surface  for  a  P  wave  travelling  in  the  unidirec¬ 
tional  fiberglass  epoxy  composite  shown  in  Fig.  1  is  an  oblate  sphe¬ 
roid.  See  Fig.  3.  Thus,  if  (x,  y,  z)  are  the  coordinates  of  a  given 

point  in  the  medium,  there  will  be  only  one  P  wave  front  passing 

•k  'k 

through  it,  corresponding  to  a  point  (S^,  Sy,  S^)  on  the  slowness 
surface  where  the  normal  is  parallel  to  the  given  direction.  There¬ 
fore,  it  is  concluded  that  for  the  displacement  components  along  the 
X,  y  and  z  axes,  u,  v  and  w,  measured  by  tracing  the  P  wave  in  the 
unidirectional  fiberglass  epoxy  composite  plate  specimen  at  point  M  in 

k  k  -k 

Fig.  8,  there  is  only  one  point  (S^,  Sy,  S^)  on  the  slowness  surface 
where  the  normal  is  parallel  to  a  given  direction  OM'  [5]. 


3 .  Directivity  Function 

The  shear  stress  Ty^  and  the  normal  stress  associated  with 

the  P  wave  reaching  the  point  M'  in  Fig.  8  are  used  to  study  their 
associated  directivity  functions.  The  asymptotic  shear  stress  and 

the  asymptotic  normal  stress  obtained  by^ substituting 

k 

Eqs .  (55)  through  (58)  into  Eq.  (1)  and  then  by  setting  equal  to 

zero  as 
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zz 


.  -(C13+C44)  *  *2  ^ 

—  IC33  [  - - -  iSy  S,  Yo 


g*2  „*2 

p  y  p 


+  S,-  +  S^"-  1)  iS^  Zo]  +  Ci3  [(-^  S 


<^44  „*2 

V 


C33  „*2 


+  —  S;‘-  1)  iS*  Yq)  +  '  '"'"p— ■— Zq]) 


-(C13+C44) 


ifc*  it 

exp{iw(Syy  +  S^z  -  t) ) 


(62) 


"y^  ~  ~2^  ‘ — ~p  y  ^  ° 


.  s«-  1)  is;z„*  (fii-  sf2 

p  y  p  ^  y 


.i|3  sf  - 1)  is:  -o)  -  is;  s:^z„, 


exp{iw(Syy  +  S*z  -  t) ) 


(63) 


The  directivity  functions  associated  with  the  normal  stress  r 
in  Eq.  (62)  and  the  shear  stress  fy^  in  Eq.  (63)  will  be  evaluated. 
Due  to  the  axial  symmetry  with  respect  to  the  zonal  axis,  the  z  axis, 
of  the  transversely  isotropic  medium,  the  values  of  the  directivity 
functions  obtained  for  the  case  Sj^  “  0  hold  for  all  values  of  Sj^. 
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Consider  the  case  of  the  point  load  acting  along  the  y  direction 
only;  that  is,  Yq  ^  0  but  Xq  -  -  0  in  Eq,  (45). 

Y 

The  directivity  function  of  the  shear  stress  Ty^  associated 
with  a  P  wave  whose  slowness  vector  is  confined  to  the  plane  x  =  0  is 
obtained  from  the  amplitude  of  the  shear  stress  in  Eq.  (63)  by  setting 
Yo  =  R  “  1; 


••^Y  ^n"C44S2  ^-Ci3  ^*2  C33  ^*2 

V  ^  ^~r'  y  ~r 


(64) 


"k 

where  is  given  in  Eq.  (57)  and  is  evaluated  at  the  point  (0,  Sy, 

k  k 

S^)  on  the  slowness  surface  for  a  P  wave,  Eq.  (14);  (0,  Sy,  S2)  is  the 

point  on  the  slowness  surface  where  the  normal  is  parallel  to  the 

given  direction  and  w  is  radian  frequency.  Similarly,  the  directivity 
Y 

function  D22  of  the  normal  r22  is,  from  Eq.  (62),  for  Yq  =  R  -  1; 


* 

nY  _  ^  C13C44  ^*2  C33C44  ^*2  ^  ^ 

^zz  ”  -  L  Sy  -  s^.  -  C]^3] 

2k  p  ■'  p 

Next,  consider  the  case  of  the  point  load  acting  along  the  z 
direction  only;  that  is,  Zq  0,  but  Xq  -  Yq  -  0  in  Eq.  (45). 


Z 

The  directivity  function  Dyg  of  the  shear  stress  Ty^  associated 
with  a  P  wave  whose  slowness  vector  is  confined  to  the  plane  x  =  0  is 
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obtained  from  the  amplitude  of  the  shear  stress  in  Eq.  (63)  by  setting 


1; 


AnwC44S*  Cii  ^*2  C33-C44-C13  ^*2 

°yz  “  — ; - ^  ^ 

^  2jr  P  ^  P 


(66) 


Similarly,  the  directivity  function  D^z  of  the  normal  stress  r^z  is, 
from  Eq.  (62),  for  Zq  -  R  -  1; 

J,Z  *n“Sz  CxiC33-Cl3C44-Cuj  s*2^i33^  S»2  ,  C33)  (67) 

2jr  p  y  p 

The  polar  diagrams  for  the  directivity  functions  of  the  shear 

Y  Y  Z  Z  ,  . 

stress  ry2  and  the  normal  stress  r2z  (^yz*  ^zz*  ^yz  ^zz  i^ 

Eqs.  (64)  through  (67))  associated  with  the  P  waves  propagating  in  the 
unidirectional  fiberglass  epoxy  composite  shown  in  Fig.  1  are  obtained 
by  substituting  Eq.  (35)  into  Eqs.  (64)  through  (67)  at  frequencies  of 
0.75,  1.50  and  2.25  MHz. 

Numerical  results  are  shown  in  Figs.  12  through  23  where  the 

angle  of  incidence  8  given  in  Eq.  (36)  is  used  to  determine  the 

*  'k 

direction  for  which  the  value  of  a  point  (0,  Sy,  S2)  on  the  slowness 
surface  where  the  normal  is  parallel  to  the  given  direction  is  thus 
obtained. 
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The  non- contact  transmitting  transducer  in  Fig.  7  is  assumed  to 
transform  an  electrical  voltage  into  a  uniform  stress;  however,  the 
non-contact  receiving  transducer  in  Fig.  7  transforms  a  displacement 
into  an  electrical  voltage.  The  approach  below  is  similar  to  that 
given  in  [ 7 ] . 

Referring  to  Fig.  7,  if  an  input  voltage  of  amplitude  V  and 
frequency  o)  is  applied  according  to 

Vi(t)  -  (68) 

the  stress  a  that  is  introduced  into  the  specimen  plate  by  the 
non- contact  transmitting  transducer  is 

a(t)  -  (69) 

where  F^Cw)  is  the  transduction  ratio  for  the  non-contact  transmitting 
transducer  in  transforming  a  voltage  to  a  stress  and  is  a  phase 
angle.  In  Eqs .  (68)  and  (69),  the  harmonic  character  of  the  signals 
is  expressed  in  the  complex  notation  where  i=»  7-1  and  only  the  real 
parts  of  these  and  subsequent  equations  should  be  considered.  Thus, 
the  amplitude  T  of  the  applied  force  is  defined  as 


T  «  Fi(o})V 


(70) 


Similarly,  if  a  stress  wave  producing  a  displacement  d  of  ampli¬ 
tude  D  and  frequency  co  that  is  detected  by  the  non-contact  receiving 
transducer  is  defined  as 

d(t)  =  (71) 

the  output  voltage  from  the  non- contact  receiving  transducer,  see 
Fig.  7,  is 

-i(c(;t  +  (^2) 

“  F2(a))De  (72) 

where  F2(c^)  is  the  transduction  ratio  for  the  non-contact  receiving 
transducer  in  transforming  a  displacement  to  a  voltage,  and  ^2  is  a 
phase  angle.  Thus,  the  amplitude  V'  of  the  output  electrical  voltage 
is 

V'  =  F2(u>)D  (73) 

The  characteristics  of  F2^(ct))  and  F2(w)  are  unknown  except  that 
the  dimensions  of  the  product  Fj^(w)F2(w)  are  [kg/m^^sec^]  . 
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5.  Steady-State  Output  Voltage  Amplitude  due  to  Multiple  Wave  Refit 
tions  in  a  Plate 


Since  the  P  wave  traced  in  the  unidirectional  fiberglass  epoxy 
composite  plate  specimen  shown  in  Fig.  8  is  travelling  in  the  y-z 
plane,  it  follows  from  Eqs .  (59),  (60)  and  (61)  that  only  the  dis¬ 
placement  components  along  the  y  and  z  axes,  v  and  w,  are  detectable 
at  the  point  M' .  Consider  first  the  point  load  acting  along  the  y 
direction  only;  that  is,  ^  0,  but  =  Zq  =  0  in  Eq.  (45). 


The  amplitude  of  the  y-component  displacement  D.^  evaluated  at 
the  point  M'  can  be  obtained  from  Eqs.  (37)  and  (60)  as 


Y  fp(Sy,S2.) 


*  *  *2  *2 
where  f^  (Sy,S2)  -  ( A^/Zir)  |  (C44//>) Sy  +  (033/^)82  -l|.  Similarly,  the 

Y 

amplitude  of  the  z-component  displacement  evaluated  at  the  point  M' 
can  be  obtained  from  Eqs.  (37)  and  (61)  as 


where  f2(Sy,S2)  -  (An/27r)  [  (013+044)/^] S„Sj 
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Next,  consider  the  point  load  acting  along  the  z-direction  only; 
that  is,  Zq  0,  but  Xq  -  Yq  -  0  in  Eq.  (45). 

Z 

The  amplitude  of  the  y- component  displacement  D.^  evaluated  at 
the  point  M'  can  be  obtained  from  Eqs .  (37)  and  (60)  as 


2  ^2  ^^y ’ 

Dv - -  Zo 

Rn 


(76) 


evaluated  at  the  point  M'  can  be  obtained  from  Eqs.  (37)  and  (61)  as 


2  f3(^y>^z) 

Dw - Z, 

Rn 


(77) 


*  *  *2  *2  . 
where  f3(Sy,  S^)  =  (Ajj/27r)  |  (C;iLl/^)^y  '^\  • 

Finally,  consider  the  point  load  acting  along  the  x  direction 
only;  that  is,  0,  but  -  Z^  -  0  in  Eq.  (45).  It  follows  from 

Eqs.  (60)  and  (61)  that  the  amplitude  of  the  y-component  displacement 
and  the  amplitude  of  the  z -component  displacement  vanish  at  the  point 
M' .  This  is  due  to  the  fact  that  the  P  wave  is  travelling  in  the  y-z 
plane . 


According  to  Eqs.  (74)  through  (77),  the  amplitude  of  a  displa¬ 
cement  component  at  the  point  M'  ,  denoted  as  ,  can  be  expressed  in 
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the  following  form  as 


if  if 

f i ( Sy , ) 

Dm< - T  (78) 

Rn 

where  T  is  the  applied  point  load  and  is  equivalent  to  either  Yq  or 

i:  if 

Zq  in  Eqs.  (74)  through  (77);  and  f^(Sy,S^)  (i=l,2  or  3)  is  determined 
from  one  of  the  Eqs.  (74)  through  (77),  depending  on  which  displace¬ 
ment  component  is  measured  and  along  which  direction  the  point  load  is 
applied. 

However,  with  the  bottom  boundary  present,  the  wave  is  reflected 
a  total  of  (2n-l)  times,  as  shown  in  Fig.  8.  Thus,  the  amplitude  of 
displacement  at  the  point  M  is  and  is  expressed  as 

2n- 1 

%  =  Qpp  (79) 

where  Qpp  is  the  amplitude  ratio  of  the  reflected  P  wave  to  the 
incident  P  wave,  as  shown  in  Fig.  5,  and  is  a  function  of  the  angle  of 
incidence  8  defined  in  Eq.  (36). 

Here,  we  have  ignored  the  effects  of  mode  conversion.  The  ampli¬ 
tude  of  the  output  voltge  from  the  non- contact  receiving  transducer  is 
V'  and  can  be  obtained  by  substituting  Eqs.  (78)  and  (79)  into  Eq. 

(73)  as 
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2n-l  ic  * 

^2^^^  Qpp  fi(Sy,S2)T 
V'  * - 

Rn 

Substitution  of  Eq.  (70)  into  Eq,  (80)  gives 


(80) 


2n-l  ^ 

F]^(cj)  F2(c^)  Qpp  ^i(^y» 

V'  -  - - - - - 

Rn 


(81) 


Introducing  the  P  wave  attenuation  constant  o  of  the  unidirectional 
fiberglass  epoxy  composite  and  a  possible  electrical  signal  amplifica¬ 
tion  factor  K,  Eq.  (81)  can  be  written  as 

2n“l  St  St  ®  ^n  /QON 

V'  -  KFi(w)F2(w)Qpp  fi(Sy,  S2)V  — — —  (82) 

Eq.  (82)  gives  the  output  voltage  amplitude  from  the  non-contact 
receiving  transducer  due  to  an  input  voltage  amplitude  V  at  the  non- 
contact  transmitting  transducer  when  the  P  wave  path  has  undergone  n 
reflections  from  the  bottom  face  of  the  unidirectional  fiberglass 
epoxy  composite  plate  specimen,  as  shown  in  Fig.  8. 
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DISCUSSION  AND  CONCLUSION 


In  the  acousto-ultrasonic  input-output  characterization  of  the 
unidirectional  fiberglass  epoxy  composite  plate  specimen,  the  angle  of 
reflection  of  the  reflected  P  wave  is  equal  to  the  angle  of  incidence 
of  the  incident  P  wave  for  each  reflection  at  either  the  top  or  the 
bottom  face  of  the  plate.  This  is  due  to  the  fact  that  the  isotropic 
plane  is  parallel  to  both  faces  of  the  plate.  However,  if  the  paral¬ 
lelism  between  the  isotropic  plane  and  the  plane  boundaries  where 
reflection  occurs  does  not  exist,  the  angle  of  reflection  is  not  equal 
to  the  angle  of  incidence.  In  this  case,  the  use  of  a  semi- inf inite 
transversely  isotropic  medium,  neglecting  the  existence  of  the  bottom 
face  of  the  plate  specimen  except  for  the  cumulative  reflections,  to 
compute  the  delay  time,  the.  displacement  and  the  directivity  functions 
becomes  inappropriate. 

It  is  observed  from  Fig,  9  that  an  increase  in  the  number  of 
reflections  n  results  in  a  minor  increase  in  the  delay  time  t^  at 
each  value  of  the  dimensionless  separation  i/h.  For  a  given  plate 
thickness  h  and  a  given  the  separation  distance  i  between  the  trans¬ 
mitting  transducer  and  the  receiving  transducer,  an  increase  in  the 
number  of  reflections  n  results  in  a  decrease  in  the  angle  of  inci¬ 
dence  dy  as  given  by  Eq.  (36),  thereby  increasing  the  travelling 
distance  R^,  from  Eq.  (37).  The  phase  velocity  of  a  P  wave  in 
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the  unidirectional  fiberglass  epoxy  composite  increases  with  decreas¬ 
ing  angle  of  8,  as  shown  in  Fig.  10.  For  the  time  delay,  the  increase 
in  the  travelling  distance  Rjj  is  partly  offset  by  the  increase  in  the 
phase  velocity  with  an  increse  in  the  number  of  reflections  n. 
This  accounts  for  the  minor  increase  in  the  delay  time  given  by  Eq. 
(38).  However,  in  an  isotropic  medium,  the  phase  velocity  of  a  stress 
wave  is  directionally  independent.  As  a  result,  the  increase  in  the 
dealy  time  tj^  caused  by  an  increase  in  the  number  of  reflections  n  is 
solely  attributed  to  the  increase  in  the  travelling  distance  in  an 

isotrpic  plate.  Therefore,  it  is  concluded  that  the  directional 
dependence  of  the  phase  velocity  of  a  stress  wave  travelling  in  a 
transversely  isotropic  medium  has  a  significant  effect  on  the  delay 
time  when  conducting  acousto-ultrasonic  testing. 

Now,  consider  the  case  of  a  given  number  of  reflections  n.  An 
increase  in  the  separation  distance  2  results  in  an  increase  in  the 
travelling  distance  R^,  as  shown  in  Fig.  8,  and  an  increase  in  the 
angle  of  incidence  8  given  by  Eq.  (36).  An  increasing  angle  of  inci¬ 
dence  8  results  in  a  decreasing  phase  velocity  C]^,  as  shown  in 
Fig.  10.  Consequently,  the  delay  time  increases  sharply  as  a 
result  of  the  increasing  separation  distance  2  for  a  given  number  of 
reflections  n,  as  shown  in  Fig.  10.  However,  as  the  separation  dis¬ 
tance  2  approaches  infinity,  the  phase  velocity  approaches  a  limit. 
Thus,  (for  i/h-*®)  the  increase  in  the  delay  time  t^^  is  due  solely  to 
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the  increase  in  the  travelling  distance  R^,  which  is  similar  to  the 
isotropic  medium  case. 

This  theoretical  investigation  provides  a  step  forward  in  the 
quantitative  understanding  of  acous to -ultrasonic  nondestructive  evalu¬ 
ation  (NDE)  parameters  such  as  stress  wave  factor  (SWF)  in  trans¬ 
versely  iostropic  media.  It  also  provides  the  potential  for  assisting 
in  the  development  of  more  efficient  and  more  revealing  NDE  schems 
utilizing  wave  propagation. 
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Laminae 


Coordinate  system  (x,y,z)  in  analysis  of  single  reflection 
problem  at  stress-free  plane  boundary  of  semi-inf inte 
transversely  isotropic  medium;  z  =  0  is  plane  boundary  where 
single  reflection  occurs,  z <  0  is  free  space,  and  x  =  0  is 
plane  of  incidence. 


••  Slowness  surface,  for  P  wave 
'Slowness  surface  for  SVwave 
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Fig.  3  Two  sheets  of  slowness  surfaces  for  P  and  SV  waves  in  unidirectional  fiberglass  epoxy 
composite  for  positive  y-z  quadrant. 


of  reflection  of  reflected  SV  wave  versus  angle  of 
ant  P  wave  in  unidirectional  fiberglass  epoxy  compo 


-component  energy  fluxof  ref  lected  P  wove 
-component  energy  fluxof  reflected  SV  wave 


Energy  fluxes  normal  to  plane  boundary,  assuming  radian  frequency  equal  to  unity, 
for  incident  P  wave  in  unidirectional  fiberglass  epoxy  composite. 


=  500 
=  300 


Fig,  13  Polar  diagram  for  directivity  function  of  shear  stress  Xyz  at  transmitting 

origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composite  at 
frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point  load  acting 
along  y-direction. 
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Fig.  14  Polar  diagram  for  directivity  function  of  shear  stress  ot  transmitting 

origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composite  at 
frequency  of  2.25  MHz  for  positive  y-z  quadrant  due  to  applied  point  load  acting 
along  y-direction. 
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Polar  diagram  for  directivity  function  Dzz  of  normal  stress  r^z  transm 
origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composit 
frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point  load 
along  y-direction. 


Fig.  19  Polar  diagram  for  directivity  function  Dy^  of  shear  stress  ^  transmitting 
origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composite  at 
frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point  load  acting 
along  z-direction . 
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Polar  diagram  for  directivity  function  of  normal  stress  T2z  transmitti 
origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composite  at 
frequency  of  0.75  MHz  for  positive  y-z  quadrant  due  to  applied  point  load  act 
along  z-d ir ec  t ion . 


Polar  diagram  for  directivity  function  of  normal  stress  transmitti 
origin  associated  with  P  wave  in  unidirectional  fiberglass  epoxy  composite  at 
frequency  of  1.50  MHz  for  positive  y-z  quadrant  due  to  applied  point  load  act 
along  Z“direc tion. 
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